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chapter 1 
INTRODOCTION 

This report describes the status of an investigation into 
the behavior of a composite structure in a fire environment. 
The investigation arises from the need to assess the damage 
a composite aircraft structure undergoes when exposed to a 
carrier deck pool fire. Towards this end, the Naval Weapons 
Center (NWC) in China Lake has undertaken a series of exper- 
iments to determine the effects of fires on structures ( 1 ). 
As the number of parameters affecting thermal behavior is 
large, and as the behavior is inherently nonlinear, an ex- 
haustive experimental study is prohibitive in cost and tine. 
To assist their experimental program, NWC has funded the 
NA7AL FOSTGRADUATT SCHOOL (NFS) to develop 5 mathematical 
model of the problem. The NWC experimental investigation, 
and the NFS analytical investigation are intended to comple- 
ment one another. 

The NFS effort seeks to quantitatively predict the extent 
of damage a composite structure undergoes when exposed to a 
fire. This objective was to be achieved through the devel- 
opment of a mathematical model of the problem. The computer 
program resulting from the mathematical model would be used 
to predict the thermal behavior of composites for a wide 
range of system parameters. The parameters character i zin g a 
system include ( 1 ) composite components (materials), ( 2 ) 
composite layup (lamina orientation), (3) composite thick- 
ness, (4) permeability, (5) pcrosity, (5) wind, and (7) fire 

intensity. The degree of validity of the model could be es- 
tablished by comparing the results of the model to experi- 

mental results. Moreover, the model would be used to sug- 
gest the most appropriate experiments. 

The overall problem of the behavior of a composite in a 
fire environment is conveniently partitioned into -wo inves- 
tigations; an investigation of thermal behavior, and an 
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investigation of how composite strength is affected by it’s 
temperature distribution. 

Mathematical models of both thermal behavior and strength 
behavior have been developed. The computer codes resulting 
from the mathematical models have been implemented on a num- 
ber of cases. Some results and brief descriptions of the 
thermal and strength analyses follow. 
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chapter 2 

THERMAL INVESTIGATION 



The thermal investigation is concerned with, the determi- 
nation of the temperature profile in a composite laminate 
due to exposure to a fire environment. 



2. 1 OBJE CTI VES 

The objectives of the thermal investigation are: 

1. development of a mathematical model for the determi- 
nation of behavior of a composite laminate in a fire 
e nvironment . 

2. to show how each of the system parameters affect the 
thermal behavior. 

3. to evaluate the extent, and rate of material consump 
tion due to combustion. 

4. to extend the Semenov combustion model to composite 
media. 

2.2 DESCRIPTION OP THE MODEL 

The principal investigator for the determination of ther 
mal response of a composite laminate in a fire environment, 
until the receipt of his Doctoral Degree in June 1932, has 
been Lcdr. C-5. Vatikiotis of the 0.3. Navy. Initaliy Lcdr 
Vatikiotis (then a Lieutenant) worked on the problem to sa- 
tisfy the research requirement towards an Engineer's Degree 
in Mechanical Engineering from NPS. The degree was granted 
in June 1980. A brief description of ^he mathematical mods 
for the thermal behavior follows. A more detailed descrip- 
tion is contained in references (2,3) . 

lie consider a composite of graphite fibers imbedded in a 
epoxy matrix. As epoxy combusts at a much, lower tsmperatur 
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-han graphite, the model assumes the epoxy has burnt off and 
all that remains is a porous graphite medium. , Figure 1 pre- 
sents a diagram of the geometry of rhe porous medium- A 
pressure differential across *he porous plane induces an air 
flow through the medium. The air flow provides two opposing 
actions; heat generation and heat transfer. The heat gener- 
ation results from the combustion which requires oxygen to 
be ever present. Heat is transferred through the porous 
medium by convective air flow. Behavior of the system de- 
pends upon which of these effects, heat generation or heat 
transfer, dominates. 

A one dimensional model of a porous medium was developed 
on the basis of the following fundamental laws: 

1. an energy balance on tha graphite, 

2. an energy balance on the air, and 

3. a mass balance on the oxygen. 

All of the heat transfer modes, conduction, convection, and 
radiation are included in the no del. In addition, the ener- 
gy balance equation for the graphite contains a heat genera- 
tion term of Arrhenius type arising from combustion of the 
graphite. Similarly, the oxygen mass balance equation con- 
tains a term accounting for rhe consumption of oxygen due to 
combustion . 

In addition to the energy and mass balance aquations, 
Darcy's law fcr the governing aquation for air flow through 
the porous graphite plate was combined with Bernoulli's 
equation to yield an equation for tha pressure differential. 
As very large changes in temperature occur during the tran- 
sient history of the system, all thermophysical properties 
are treated as temperature (and hence time) dependent. The 
combustion model was patterned after the Semenov model, and 
the reaction was taken as producing only carbon dioxide. 

The mathematical model of the thermal problem consists of 
three nonlinear, coupled, partial differential equations in 
time and space. These field equations were transformed, via 
a Galerkin formulation of tha finite element method, into a 
set of 3 xn nonlinear, coupled ordinary differential equa- 
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Figure 1. Geometry of the Model 



wions in time, where r. is -he number of spatial nodes in the 
finite element discretization. These equations, upon numer- 
ical integration by the FranJce (4) modified Gear method for 
stiff systems, gives the graphite nemperature, the air temp- 
erature, and the oxygen concentration at each nodal point as 
a function of time. In addition to these response varia- 
bles, a large number of other system parameters are calcu- 
lated. These include air density, porosity, permeability, 
Reynold's number, pore velcciry, and pressure. In this way 
a full characterization of the system is known during the 
transient history of the problem. 

2,3 GE52aiLI?ATI^ OP THE MODEL 

On the basis of his NPS research activity, Lcdr. 7ati- 
kiotis was accepted to the NPS Doctoral program. He was 
permitted to conduct the research required for the Doctoral 
Degree at his next duty station, the David Taylor Naval Ship 
Research and Development Center (NSRDC) . During the period 
of June 1980 to June 1982, Lcdr. Vatikiotis made significant 
improvements to the thermal model, and performed a number of 
computer runs to determine the affects of design and envi- 
ronmental parameters on system behavior. A brief descrip- 
tion of his work is presented here. For a more detailed ac- 
count of the work, the reader is referred to the Doctoral 
dissertation (5) . 



2.3.1 Modul ariz ation of the Computer Progr am 

The computer program generated in the Engineer's Degree 
research effort was put into modular form in order to more 
easily accommodate any changes to the model. At present, 
the computer code consists of 25 subroutines. 

2.3.2 Tra nsient Porosity and Pe rmeability 

In the earlier model, porosity and permeability were 
treated as uniform through the medium, and constant with re- 
spect to time. The model was updated to account for an ini- 
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•cial r.cn uniform disrriburion of graphite filaments. More- 
over as graphite is consumed by combustion, both the porosi- 
ty distribuxion and permeability distribution vary with time 
and position. The model was modified to account for the 
changing air density, and porosity by invoking the continui- 
ty equation. Combining the continuity equation with Darcy's 
law resulted in a fourth nonlinear second order field equa- 
tion. This equation was solved numerically by the "shooting 
me xhcd". 



2.3.3 Combu st io n 

In the earlier model the product of the graphite combus- 
tion was assumed to be carbon dioxide only (i.s., a simple 
reaction). The literature shows that carbon monoxide be- 
comes a byproduct, along with carbon dioxide, at higher com- 
bustion temperatures. An empirical relation for the propor- 
tion of CO and C02 as a function of temperature was 
incorporated into the model. 



2,3.4 Arrheni us Sx pre ssion 

As Frank-Kamenetskii (6) suggested that carbon combustion 
is net a first order reaction, the model was generalized to 
accommodate fractional order reactions. This modification 
was accomplished by raising the oxygen concentration term in 
the reaction rate expression to a fractional power. A ser- 
ies of analyses were performed with fractional exponents 1/2, 
2/3, and 3/4. The analytical results changed markedly with 
exponent. Hence, the accuracy of any analytical model de- 
pends on the accuracy of the experimentally derived exponent. 



2,3.5 Transient Plate Thic kn ess 

Accounting for the change in porosity due to graphite 
consumption also permitted the determination of the laminate 
thickness with time. 
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2.3.6 Porous Model 

Oriqinally, the porous medium was modeled as cylindrical 
fibers. The model has been modified so than it can treat 
either spherical particles or cylindrical fibers. This op- 
tion permits the user to consider the problem of combustion 
of graphite particles. 



2.3.7 Heat Transfer Co ef fic ient 

Alternate models for internal heat transfer coefficient, 
surface area per unit volume, and average pore diameter are 
available for selection cf the investigator. in addition, 
the pore diameter is based on the hydraulic diameter. 

2.3.8 Surface R ec ession Pr o blea 

In the case of combustion, the reaction invariably moves 
to the source of the oxygen, namely the air inlet surface. 
This results from the fact that in a relatively short time, 
all cxygen in the interior is consumed. Dnce all of the in- 
terior oxygen has been consumed, all subsequent combustion 
takes place at the air inlet surface where there is always a 
continuous supply of cxygen. At this stage in the transient 
history, the problem degenerates into a "surface recession" 
pr cblem. Since there is no oxygen except at the inlet sur- 
face, the oxygen balance equations must be removed from the 
analysis, otherwise severe numerical difficulties arise. 

The removal cf these equations from the analysis is accom- 
plished automatically in the computer program, so that the 
user need not be concerned with the transformation to a sur- 
face recession problem. 



2.3.9 Initial C on d ition s 

The arbitrariness of selecting initial conditions was 
cverccme by incorporating into the analysis, a subroutine 
which generates a transient thermal environment. If the 
user selects this mode for the initial conditions, a heat 
flux and it's duration time are specified by the user. 
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2.3.10 Bou ndary Con diti ons 

The model was given more flexibility for applica-icn to a 
wider class of problems by allowing a choice of boundary 
conditions. The user may select any of the following three 
boundary conditions: 

1. insulated surfaces and Dankwert's B.C.'s for the air 
temperature and oxygen concentration. 

2. radiarion at surfaces and Dankwert's B.C.’s for the 
air temperature and oxygen concentration. 

3. convective and radiative heat transf'='r at both sur- 
faces, air temperature at inlet surface equal to am- 
bient temperature, oxygen concentration at -he inler 
surface equal to ambient concentration, and Dank- 
wert's B.C.'s at the exit surfaces. 

2.3.11 numeric al In st abilit ies 

During certain times in the transient history, numerical 
ins t abil it ies were encountered. This is particularly true 
when the oxygen concentration is going to zero in the inte- 
rior of the medium. The numerical instability appears in 
the results as an oxygen concentration oscillating about 
zero. A subroutine was added to the program which dampens 
the oscillation significantly. 



2.3.12 Impr oved Accur ac y 

In the original formulation of the problem, the nonlinear 
heat generation term was treated as a "forcing" term. A 
significant reduction in computational effort resulted when 
this nonlinear term was treated as a part of the field oper- 
ator. This improvement in computational effort is attribut- 
ed to the fact that as an operator, the generation term en- 
ters the analysis as part of the Jacobian matrix of the 
system. As the Jacobian matrix plays a prominent role in 
the numerical integration, it is no wonder that a signifi- 
cant improvement occurs. 
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2.3.13 O pdatinq Th eraophys i c al Progerties 

In the original formulation, all temperature dependent 
properties were updated outside of the time integration 
loop, at the arbitrary discretion of the aser. The formula- 
tion was modified so that all temperature dependent proper- 
ties are updated continuously and automatically within the 
integration loop. The result is improved accuracy and a re- 
duction of computational effort. 

2.3.14 C on d uc tivity 

The expression for thermal conductivity of the porous 
medium was changed to one that extended the range of porosi- 
ties which could be treated. The magnitudes of the internal 
parameters, such as pore velocity, etc., are checked to ver- 
ify that the empirical expressions used to calculate these 
parameters remain valid. 

All of the modifications briefly described above had the 
effect of (1) extending the modal to a more general stats, 

(2) provide the analysis with more accuracy, and (3) reduce 
the computational effort required in the solution of a prob- 
lem. In the next section some of the results obtained by 
Lcdr. yatikiotis are briefly described. 

2.4 DISS ERT ATION RESULTS 

A number of computer analyses were performed in order to ob- 
tain results which would provide some understanding of the 
effects of each of the system parameters on system behavior. 
In order to establish the validity of the mathematical mod- 
el, a comparison of results frcm computer analyses with ex- 
perimental results was made whenever possible. In such 
cases there was fair agreement between analytical and 
experimental res'ilts. There ace three reasons why excel- 
lent agreement between experimental and analytical results 
should not be expected. First, combustion experiments by 
their very nature are difficult to perform and control. 
Thus there is a large uncertainty in the reported results. 
The relatively small number of combustion experiments at- 
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tssxs to t.heir difficulty. Secondly, zhere is small likeli- 
hood that the values of many system paraneters is known to 
any high degree of certainty. Thus, for example, the expo- 
nent on oxygen concentration in the Arrhenius expression has 
not been accurately determined; and yet a small change in 
the exponent yields significant changes in the quantitative 
results. Thirdly, there is the inherent uncertainty associ- 
ated with every complex mathematical model. Any number of 
simplifying assumptions are made; and of course there are 
choices in the selection of competing relations for calcu- 
lating properties. All of these considerations suggest that 
the results be regarded as descriptive rather than precise. 
As ncne of the analysis results contradicted experimental 
results, there is confidence that the mathematical model is 
a good one. Figure 2 presents a comparison of the mathemat- 
ical model and the experimental results of Kolodstev (7) . 
The data for the analysis is given in Table 1. 

A series of computer runs were performed to isolate the 
effects of a single system parameter on system behavior. A 
brief description of each of the studies follows. In each 
case, the boundary conditions with insulated surfaces was 
used. The numerical results given in this report should be 
considered approximate rather than precise, 

2.4.1 Ef fec ts of Fire 5n vironm ent 

In this study, the effect of fire intensity on a compos- 
ite plate was investigated by subjecting a plate to two 
heat flux conditions. The first heat flux condition results 
in combustion, and the second heat flux ccndition results in 
ext ingui shment. Table 2 presents the data describing the 
composite of this study. 

When a heat flux of 30, 000 3tu/ft-hr is applied to the 
inlet surface for 15 seconds, and is then removed, the plate 
undergoes combustion. Figure 3 shows the graphite and air 
temperatures as functions of space at selected times. Fig- 
ure 4 plots oxygen concentration versus space at selected 
times, and Figure 5 shows the thickness of the plate as a 
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TABLE 1 Geometry o-f porous medium and ambient conditions 

-For simulating Kolodstev's C 7 1 experiment. 



• Particle Shape: 


spherical 


• Particle Diameter ( d ) : 


. 126 in 


Unit Cell Thickness 


. 126 in 


• Spatial Thickness of Porous Medium ( L. ) : 


7.5 in 


Porosi ty < p ) 


. 476 


Permeability (m> 


.103x10^ ft^ 


Bulk Thermal Conductivity of Carbon (k^); 


86.0 Btu/ft-hr-F 


Bulk Specific Heat of Carbon ( ) : 


.231 Btu/lbm-F 


Bulk Density of Carbon (p ): 


70.3 Ibrn/ft"^ 


Thermal Emissivity of Particles ( E ); 


. 9 


Ambient Temperature ( T*) ! 


80.0 deg-F 


• Ambient Pressure (f^): 


14.7 psi 


• Pressure Differential <aP) ; 


-.29 psi 


• Ambient Oxygen Concentration <C^») : 


.0172 Ibm/ft^ 



•Known parameter for Kolodstev's experiment. 
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function of time. Figure 4 shows rhar the graphite and air 
temperatures differ only at the air inlet surface (X/L=0.) 
until about t=360 sec. Thereafter the air temperature is 
significantly less than the graphite temperature for a con- 
siderable range of X/L. This can be misleading because the 
thickness L has changed signif ic anrly during the evolution. 
Figure 5 shews that the plate thickness at 517 seconds is 
only 10 percent of it's original thickness. Thus the pene- 
tration of "cooler" air is .0025 inches. At 20 seconds, the 
plate thickness is .25 inches and the cooler air penetrates 
the plate to X/L=0.1, or .025 inches. Thus when looking at 
Figure 3, one must keep in mind that L, the rhickness of the 
plate, is changing due to combustion. 

An important observation which can be made for ■^his com- 
bustion case cones from Figure 4. Figure 4 shows that the 
oxygen in the interior of the plate is very quickly con- 
sumed. All cases which have gone ro combustion have shown 
this trait, the interior oxygen is quickly depleted. At 
time equal ro 20 seconds, the problem is about to become a 
"surface recession" problem, in which all reacrion takes 
place at the air inlet surface. Figure 5 shows that the 
rate of material consumption bee ernes ever increasingly 



greater with time after 


surface recession 


begins. For 


“this 


particular case we note 


that the temperatu 


re increased 


from 


a maximum of 900 degrees 


Fahrenheit when t 


he heat flux 


was 



removed ro a temperature slightly less than 3000 degrees 
Fahrenheit at 51t seconds. 

When a heat flux of 30, 000 3tu/ft-hr is applied to the 
inlet surface for 14 seconds and then removed, the plate 
goes to extinction in less than a minute. Figure 6 shows 
the manner in which the temperature distribution goes to the 
uniform ambient temperature. Se note that -^he peak tempera- 
ture while continuously decreasing in magnitude with time, 
moves from the inlet surface to the exit surface. Figure 7 
shows that as the plate cools, the oxygen concentration in- 
creases until it reaches the uniform ambient condition. 
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TABLE 2 Geometry o-f porous medium and ambient conditions 



■for the example problems in 


Section 2.4.1 



Particle Shape: 


spherical 


Particle Diameter (d): 


. 005 i n 


Unit Cell Thickness ( 0)J 


.005 in 


Spatial Thickness o-f Porous Medium ( L > s 


.25 in 


Porosi ty < p ) 


. 476 


Permeability ( nn ) 


-q 2 

.162x10^ ft^ 


Bulk Thermal Conductivity o-f Carbon ( k^) : 


86.0 Btu/ft-hr-F 


Bulk Speci-fic Heat o-f Carbon ( ) : 


.231 Btu/lbm-F 


Bulk Density o-f Carbon < > s 


70.3 lbm/-ft^ 


Thermal Emissivity o-f Particles ( € ): 


.9 


Ambient Temperature (T ): 


30.0 deg-F 


Ambient Pressure ( P^) : 


14.7 psi 


Pressure Di -f -f erenti al (Ap) : 


-.35 psi 


Ambient Oxygen Concentration 


.0172 Ibm/ft^ 
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FIGURE 3 Sustained combustion resulting from a heat 

flux at the x = 0 surface. 
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FIGURE 4 Oxygen response during sustained combustion. 
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FIGURE 5 Change in thickness of porous medium during 

surface recession phase. 
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FIGURE 6 Response of temperature when heat flux at x 

insufficient to produce sustained combustion 
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FIGURE 7 Response of oxygen as porous medium cools 

to ambient temperature . 



The point to be made here is that a slight change in rhe 
duration of hear flux can determine whether damaging combus- 
tion or safe exrinction occurs. Here we have seen that 
30,000 3tu/ft-hr for 15 seconds resulrs in a combustion 
which completely consumed the plate. The same plate goes *o 
extinction if the heat flux is terminated after 14 seconds. 

2.4.2 Eff ec ts of Perme abil ity 

Permeability, a measure of flew conductance, is one of 
the composite design parameters which most affects thermal 
behavior. In order to isolate the effect of permeability on 
behavior, the geometry of the filament diameter and the 
width of the unit cell were varied in such a way as to main- 
-ain a constant porosity for all analyses in this study. In 
this way, the quantity of air around a filament at any in- 
stant is fixed, but the rate of air passing a filament in- 
creases with increasing permeability. Table 3 lists the 
data for a number of analyses which were used in the perme- 
ability study. Figure 3 presents tempera*: ure versus perme- 
ability for three values of pressure differential. The 
temperature hers is the temperature above which the plate 
will experience combustion, and below which the plate will 
extinguish. The "ignition" temperature for each case is ob- 
tained by trial and error iteration on a uniform initial 
temperature. After a number of runs, we obtain two nearly 
equal uniform initial temperatures, one leading to combus- 
tion and the other extinguishment. Their average is taken 
as the "ignition" temperature. 

Figure 8 shows that for a giver, permeability, the temper- 
ature at which combustion occurs increases with increasing 
pressure differential. An increase in the pressure iiffer- 
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an increase in the ignition temperature. We also note that 
the rate of change of ignition temperature decreases with 
increasing permeability. 
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TABLE 3 Geometry of porous medium and 

for the permeability analysis 


ambient conditions 
in Section 2.4.2 



Particle Shape: 


spherical 


Particle Diameter ( d ) : 


(various) 


Unit Cell Thickness 


(various) 


Spatial Thickness of Porous Medium ( L 


1.0 in 


Porosi ty < p ) 


. 476 


Permeability m ) 


(various) 


Bulk Thermal Conductivity of Carbon < : 


86.0 Btu/ft-hr-F 


Bulk Specific Heat of Carbon ( d: ) : 


.231 Btu/lbm-F 


Bulk Density of Carbon > : 


70.3 Ibm/ft'^ 


Thermal Emissivity of Particles (€ >: 


. 9 


Ambient Temperature (1^): 


80.0 deg-F 


Ambient Pressure ( P): 

CO 


14.7 psi 


Pressure Differential s 


( various) 


Ambient Oxygen Concentration (c£^) : 


.0172 Ibm/ft^ 
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2.4.3 Ef fe cts of P oros ity 

The affect of porosity, the a 
unit volume of space, on thermal 
9. The data for the porosity st 
Figure 9 shows ~hat for a fixed 
perature increases with increasi 
shows rhat increasing the air fl 
fixed temperature, or stating it 
creasing the air flow through a 
to combust. 

However, given a fixed permea 
ature decreases with increasing 
this is that increasing the poro 
parameters results in decreasing 
um. We note that the rate of ch 
diminishes with porosity. That 
temperature with porosity become 
sit y . 

2.4.4 Ef fe c ts of Thickne ss 

Figure 10 shows the effect of 

behavior. Table 5 is the data f 
most imporzant feature of Figure 
crease in ignition zemperature w 
for the smaller thicknesses from 
An explaination for this behavio 
paramezers except thickness are 
thickness while keeping the pras 
rssulzs in decreasing the pressu 
in pressure gradient causes a ie 
of the air through the plate. T 
results in a decrease of air flo 
veczive heat transfer diminishes 
ni~y for blowout. Hence an incr 
suit in a lower ignition tempera 

Anozher explaination for deer 
with increasinc thickness is as 



easure of void volume per 
behavior is shown in Figure 
udy is given in Table 4. 
porosity, the ignition tem- 
ng permeability. Again this 
ow will cause blowout at a 
a little differently, in- 
plate makes it less likely 

bilizy, the ignition t-mper- 
porosity. The reason for 
sizy while fixing all other 
air flow through the medi- 
ange of ignition temperazure 
is, zhe decrease of ignition 
s less wi‘-h increasing poro- 



plate thickness on thermal 
or the thickness szudy. The 
10 is the very large de- 
ith thickness, especially 
0.25 inches to 2.0 inches, 
r is the following. If all 
fixed, then increasing the 
sure differential constant 
re gradient. The decrease 
crease in the pore velocity 
hus, increasing thickness 
w. With less air flew, con- 
and there is l^ss opportu- 
ease in thickness should re- 
ture . 

easing ignizion temperazure 
follows. The above sentence 
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TABLE 4 Geometry o-f porous medium and ambient conditions 



for the porosity analysis 


in Section 2.4.3 



Particle Shape: 


spherical 


Particle Diameter ( d ) : 


(various) 


Unit Cell Thickness 


( vari ous) 


Spatial Thickness of Porous Medium (L)s 


1 . 0 i n 


Porosi ty ( p ) 


(vari ous) 


Permeability (m) 


(various) 



Bulk Thermal Conductivity o-f Carbon ( : 36.0 Btu/ft-hr-F 



Bulk Specific Heat of Carbon (C^): 


.231 Btu/lbm-F 


Bulk Density of Carbon (p,.^* 


70.3 Ibm/ft^ 


Thermal Emissivity of Particles (€ )s 


.9 


Ambient Temperature ( T^) : 


30.0 deg-F 


Ambient Pressure ( ? ) ; 


14.7 psi 


Pressure Differential (aP) : 


-.5 psi 


Ambient Oxygen Concentration (<p : 


.0172 Ibm/ft^ 



25 




26 



FIGURE 9 --- Porosity POROSITY 

analysis. 



TABLE 5 — - Oeometry o-f porous medium and 

for the thickness analysis in 

Particle Shape: 

Particle Diameter (d ): 

Unit Cell Thickness 

Spatial Thickness of Porous Medium ( L ): 
Porosi ty < p > 

Permeability (m) 

Bulk Thermal Conductivity of Carbon ( k^) : 
Bulk Specific Heat of Carbon ( ) s 

Bulk Density of Carbon ( Pg ^ s 
Thermal Emissivity of Particles (€ ): 
Ambient Temperature ( T^) : 

Ambient Pressure ( f^) : 

Pressure Differential (Ap> : 

Ambient Oxygen Concentrati on ( ) : 



ambient conditions 
Section 2.4.4 

spheri cal 
(various) 
(various) 
(various) 

.476 

(various) 

86.0 Btu/ft-hr-F 
.231 Btu/lbm-F 
70.3 Ibm/ft^ 

.9 

ao.O deg-F 
14.7 psi 
-. 5 psi 
.0173 Ibm/ft^ 
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POROUS MEDIUM THICKNESS ( in ) 



can te turned arcund to ’’decreasing the thickness will re- 
sult in increasing the ignition temperature". Now the ~hin 
r.er the plate, the less distance the air traverses in pass- 
ing through the plate. Since the air temperature increases 
as it passes through the plate, thin plates have a larger 
fraction of their region with cool air than do thick places 
As thin plates are more effectively cooled by the air flow, 
higher ignition temperatures are needed for combustion. 

The three curves again show than increasing the pore ve- 
locity by increasing the permeability, leads to blowou-, an 
hence greater air flow (or, equivalently, permeability) 
leads to higher ignition temperatures. 

Another striking feature of Figure 10 is the rapid de- 
crease in the rate of change of ignirion temperature with 
chickness. A large change in the ignition -emperatura can 
be effected by a small change in the thickness when the 
plate is thin. On rhe other hand, very little change in th 
ignition temperature can be effected if the plate is very 
thick. 

2. U . 5 Ef fec ts of Pressur e D if f e rent ial 

Each of the previous studies sought to determine the ef- 
fect of a design parameter on system behavior. Thesa param 
eters are, generally, under the control of the structural 
engineer, and thus can be selected to minimize combustion. 
In contrast, there are external parameters such as pressure 
differential (possibly arising from wind over one surface) 
which are difficult, if at all possible, to control. 

Table 6 presents the data for the pressure differential 
study, and Figure 11 presents the results of the study. Th 
results show that the ignition tenpe nature increases with 
increasing pressure differential, and that the rate of in- 
crease decreases with pressure differential. As in the pre 
vious studies, it is again evident that an increase in 
permeability results in higher ignition temperatures. The 
three curves in Figure 11 show that the effect of pressure 
differential is much greater for larger values of permeabil 
ity. 
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2.4.6 Effects of B oa ndar y C on di tio ns 

All of the the previous parameter studies wers performed 
using insulated boundary conditions at the inlet and exit 
surfaces. In this study, the effect of hear transfer at the 
boundaries was investigated. Figure 12 compares the behav- 
ior of plates with insulated boundaries (curve 3) to plates 
with heat transfer at the boundaries (curve A). Although 
the ignirion temperatures for the plates wirh heat transfer 
are always greater than the ignition temperatures for the 
insulated plates, two regions are distinguishable. In the 
region with permeability greater than 0.2x10-8 in. , the ig- 
nition temperatures for insulated plates are but slightly 
less -han for the heat transfer plates. In the region with 
permeability less than 0.2x10-8 in. , there is a far mors 
pronounced difference between insulared plates and plates 
with heat transfer. The results show that as permeability 
decreases, the ignition temperature decreases in bonh cases, 
but the rate of decrease is larger for the insulated case, 
and thus the difference between ignition temperatures be- 
comes increasingly greater. Finally, at some small value of 
permeability, the ignition temperature for plates with sur- 
face heat transfer reverses and begins increasing rapidly 
with decreasing permeability. 

An increase in the ignition temperature with decreasing 
air flow (permeability) for plates with heat transfer bound- 
ary conditions could be explained if the heat generation by 
the exothermal reaction , and not the heat transfer mechan- 
ism, was the dominating mechanism. In this case the reduc- 
tion of oxygen resulting from a decrease in permeability 
would extinguish the combustion by ’’chok-ing”. In the chok- 
ing region, ignition temperatures should increase with de- 
creasing permeability. 

Apparently, porous plates with heat transfer boundary 
conditions have a minimum ignition temperature at a particu- 
lar value of air flow, say U. Increasing the air flow above 
U results in combustion extinguishment by blowout and hence 
greater igni-^ion temperatures. This is the case when heat 



30 



TABLE 6 



Geometry o-f porous medium and 

■far the pressure di f -f erenti al 
Section 2.4.5 

Particle Shape: 

Particle Diameter ( d ) : 

Unit Cell Thickness 

Spatial Thickness o-f Porous Medium < L ) s 
Porosi ty < p ) 

Permeability (m) 

Bulk Thermal Conductivity o-f Carbon ( k^) : 
Bulk Sped -fie Heat o-f Carbon ( C,. > : 

Bulk Density o-f Carbon t 

Thermal Emissivity o-f Particles (€ ); 
Ambient Temperature < J 
Ambient Pressure < P ) : 

Pressure Di -f -f erenti al (Ap) : 

Ambient Oxygen Concentration (Cb ) ; 



ambient conditions 
analysis in 

spherical 

(various) 

(various) 

1.0 in 
. 476 

( various) 

86.0 Btu/-ft-hr-F 
.231 Btu/lbm-F 

70.3 lbm/-ft^ 

. 9 

80.0 deg-F 
14.7 psi 

(various) 

.0172 Ibm/ft^ 
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PERMEABILITY Iff 



transfer is the dominant mechanism. On the other hand, de- 
creasing the air flow below (J results in combustion extin- 
guishment by choking and hence greater ignition tempera- 
tures. This is the case when the chemical reaction is -he 
dominant mechanism. 

2.4.7 Effects of Pore Velo city 

Consideration of the results of the previous studies on 
rha effects of individual parameters on system behavior sug- 
gesrs that the most important parameter affecting behavior 
is -he pore velocity. With this in mind, the results of all 
computer analyses (47 problems) are shown in Figure 13. The 
curve through these data points is a fair fit, and thus 
there is reason to believe that the ignition temperature may 
be obtained, with fair accuracy, by the equation of the 
curve. 

This concludes the summary of the work done by Lcdr. 7a- 
tikictis at WSRDC for his Doctor of Engineering Degree. The 
dissertation, a more general treatment of the combustion 
problem, presents a more detailed discussion of the model, 
numerical aspects of the formulation, the results, and con- 
clusions. In addition, there is a discussion of a general- 
ized Semenov model for combustion of porous media, as well 
as an extensive list of references. 



2.5 NPS ACT IVITY 

During the past year, about half of the funded research 
activity at NPS was devoted to the thermal investigation, 
and the other half was given to the strength investigation. 
In this section, a summary of the effort of the thermal in- 
vestigation at NPS will be described. 



2.5.1 Resta rt Capability 

One of the input parameters 
gram (CARBGT) is the time that 
ysis to stop. It is difficult 
for this parameter. Selecting 



n the thermal 


analysis 


pro- 


he user would 


like t h e 


dfi dl 



to select a reasonable value 
too large a time is wasteful 



34 



0 

9 

I 




L 



U 

q 

UJ 

> 

cn 

UJ 



35 



FIGURE 13 "Minimum initial temperature" as a function of pore velocity. 



of ccmputer effort, while roo small a tine could require re- 
doing the problem. To remedy this situation, -a subroutine 
was added to the program to permit restart. The restart 
feature permits the program to start from the results of a 
previous analysis. All of the input data necessary for the 
restart analysis is obtained from the final output of the 
previous analysis. Restart can be invoked in the analysis 
of a problem as many times as desirable. In this way, a 
complere analysis of a problem may consist of several time 
contiguous analyses. 



2.5.2 Ef fec ts of I n it ial Conditions 

There are two ways to input initial conditions. A way to 
obtain reasonable initial conditions on temperature is to 
specify a heat flux and a duration time, as described in 
section 2.3.9. On the other hand, if the user is interested 
in the effects of the initial temperature on thermal behav- 
ior, he may input any initial temperature distribution. 

This was the method used in this investigation. Here we ob- 
rained the effects of linear initial temperature distribu- 
■^ions on ignition temperature in the following way. The 
slope of the linear temperature and the average temperature 
was specified. The same problem was run as many times as 
needed to obtain the average temperature below which extinc- 
tion occurs and above which combustion occurs. Figure 14 
presents the results of this study. Two features of this 
fiaure are (1) the maximum ignition •temperature occurs for a 
negative temperature slope, and (2) the ignition tempera- 
tures are greater for the negative temperature slopes. Fig- 
ure 14 shows that combustion is less likely to occur for 
plates with the higher temperature on the X/L=0. surface. 
This is a reasonable expectation since the cooling air is 
most effective at it's entry into the medium. 
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temperature difference IF) 



2.5.3 Gra phic s C apab ility 

Each computer analysis produces a vast quantity of numer 
ical output. Thus, at preselected instances of time and ar 
each spatial nodal point, the program reports: (a) graphite 

temperature, (b) air temperature, (c) oxygen concentration, 
(d) reaction rare, (e) filament diameter (f) porosity, (g) 
permeability, (h) pressure, (i) pressure gradient (j) pore 
velocity (k) Reynolds Number, (1) convection coefficient, 
and (m) area to volume ratio. In addition, plate thickness 
is reported at each time. This data must be assimilated by 
-he program user if full benefit of the analysis is to be 
obtained. As a typical analysis consists of 30 to 60 time 
steps and 25 no 40 spatial nodal points, a typical output 
contains anywhere from 10, 000 tO 30,000 numbers. It is dif 
ficult, if a- all possible, to assimilate this quantity of 
numerical data, and therefore graphical output is more than 
just desirable, it is necessary. To this end, araphics ca- 
pability was incorporated into the program. A variety of 
graphics options are now available. The user may select an 
combination of the following output modes: 

1. numerical results of all system parameters, 

2. curves of graphite temperature and oxygen concentra- 
tion as functions of space at preselected time in- 
stances, 

3. three dimensional surfaces of graphite temperature, 
oxygen concentration, and reaction rate as functions 
of time and space, 

4. contour graphs of graphite temperature, oxygen con- 
centration, and reaction rats. 

Any combination of these options is available a-^. a CRT ter 
minal or as hard copy output. 

Two programs exist for the surface and contour graphics. 
GRAF3D produces surfaces and/cr contours from the actual 
numerical output of the analysis program, while GRAF3E 
achieves smoother graphics by passing the outpu-^ from CAH3G 
through a cubic spline int erpola ticn subroutine (IMSL sub- 
routine IBCIEU) . 



from GRA?3D or 



After inspecting the graphical output 
GRAF3E, it may be desirable to have a mors detailed surfac 
or ccntour of a local subregion of the time-space domain. 
This may be done by recalling either graphics program, and 
specifying the upper and lower limits of the rirae and spac 
subregion of interest. Figures 15 through 20 are sample 
output of surfaces and contours from GRAF3D. Figures 21 
through 26 , are samples of subregion graphics of the same 
analysis as Figures 15 through 20. 

This concludes the discussion of the research on the 
thermal investigation. A description of the strength inve 
nigation follows. 
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Figure 15. Temperature Surface from GFuAFSB 
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Figure 16. Concentration Surface from GRAF3D 
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Figure 17. Reaction Rate Surface from GRAF3D 
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Figure IS . Temperature Contours from GRAF3D 
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igure 19. Concentration Contours from GRAF3D 
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Figure 20. Reaction Rate Contours from GRAF3D 
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Figure 21. Subregion Surface of Figure 15 
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Figure 22. Subregion surface of Figure 16 
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Figure 23. Subregion Surface of Figure 17 
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Figure 24. Subregion Contour of Figure 18 
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Figure 25. Subregion Contour of Figure 19 
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Figure 26, Subregion Contour of Figure 20 



chapter 3 

STRENGTH INVESTIGATION 

When a composite laminate is subjected to a thermal dis- 
tribution, thermal stresses will arise if thermal expansion 
is restrained. The stresses in a laminate are the result of 
-hermal and load stresses. The strength of a laminate de- 
pends upon the stress state and the failure criterion go- 
verning failure. 



3 . 1 OBJE CTI VES 

The objectives of the strength investigation are; 

1. to determine the stresses in a laminate due to any 
combination of inplane loads, bending loads, and 
"thermal loads", 

2. to determine the strength of a laminate resulting 
from these stresses. 

3, 2 description OF THE M^BL 

A detailed formulation and description of the mathemati- 
cal model is presented in a technical report (8). A brief 
summary follows. 

The model is restricted to symmetric laminates. This is 
hardly a limitation on the modal because 

1. nonsymmetric laminates are extremely uncommon, and 

2. symmetric geometry does not restriot the model to 
symmetric temperature distributions. 

The formulation begins with the formation of the stiff- 
ness matrix of a lamina in the natural coordina'^es of the 
lamina. The stiffness matrix of a laminate is obtained from 
the lamina stiffness matrices by coordinate transformations 
from the natural axes of a lamina to the laminate coordinate 
system. These are the laminate stress-strain relations. 

The strain-displacement relations {Kir ohhof f-Love theory) 
of thin plate theory are combined with the stress-str ain re- 
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laticns zo give zhe lamina-e force-displacement relations. 
The stresses in each lamina can be obtained from the forces 
through the force-stress relations. Temperature effects en 
ter into the analysis through the stress-strain relations, 
which in turn produced "thermal load" terms in the force- 
stress relations. 

Once the lamina stresses are known, the task is to deter 
mine which, if any, of the lamina have failed. The general 
ized Tsai-Wu criterion is used as the failure criterion for 
the lamina. The failure of one or more lamina does not im- 
ply that the laminate has failed. The snrength analysis 
de-ermines the ply by ply sequence of failure. After a ply 
(lamina) has failed, the ply is removed, and the stress- 
strain relations for the modified laminate are generated. 
This process is repeated until all lamina have failed. If 
all cf the plies have failed, the laminate has failed. If 
the loading and temperature are such that all of the lamina 
have not failed, then the laminate has not failed. In this 
case the task is to determine "how far" the laminate is fro 
failure. This is obtained from calculating the "load fac- 
tor". 

3.3 HODIFICATIONS OP THE MODEL 

The original formulation of the model was restricted to 
symmetric temperature field problems. In addition, the 
analysis was restricted to problems for which only the load 
factor could be determined. These restrictions have been 
eliminated in the present model cf the strength analysis 
program (STRSN) . 

3.3.1 Nonsymm etr ic Tem perat ure Fields 

The strength analysis program (3T3EN) was modified to 
handle nonsymmetric temperature distr ibu"^ions which result, 
for example, when one surface cf the laminate is exposed to 
a fire. Some of the equations presented in reference (8) 
are changed as fellows; 
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3.3.2 Ad di ti on al F a il ur e Fa ctors 

In it's earlier form, STHEM provided the load factor to 
failure for a specified temperature distribution. Tbe pro- 
gram was modified to handle "temperature factor", and "com- 
bined factor" problems as well. each of these problems wil 
be briefly described. 

In the load factor problem, the magnitude of temperature 
in each ply is specified explicitly, and the relative magni 
nudes of the laminate forces and bending moments are given. 
If the laminate load vector (Mx, Ny,Nxy,dx, My, Mxy) is speci- 
fied as ( 1 , , . 5, 0. , 0. , 1 0. ,0 ,) , and the load factor is derer- 
mined zo be R, then for the given temperature field, ■'•he la 
minate will fail when Nx = 1.x?., Ny=0. 5xR, Sxy=0.xR, Mx = 0.xR, 
My = 10.xR, and Mxy=0.xR. Reference (8) shows that rhe load 
factor appears as -he scalar multiplier of the "load stress 

In the temperature factor problem, •rhe laminate load mag 
uirudes are specified exactly, and only the relative magni- 
tudes of the ply temperatures are known. In this case, the 
problem is to determine the factor which -when multiplied by 
the relative lamina temperatures gives the temperature fiel 
at laminate failure. 

Finally, we consider the combined factor problem. Hers 
the relative maanitudes of both the loads, and ply tempera- 
tures are specified, and the task is to determine the "com- 
bined factor" which, upon multiplication with each or the 
relative loads and temperatures, gives the loads and temper 
azures at which the laminate will fail. 

In each case, the Tsai-Wu failure criterion reduces to 
the guadratic in the factor R, 

aR^ + 3R •*• Y = 0 

•The coefficients a, 3 / and Y depend upon which R factor is 
being determined. In the load factor case, a ,3 , and Y / ar 
the a, b, and c of Equations (55) through (67) of reference 
(8) , respectively. In the temperature factor case, a is b, 
is a, and y is c. In the combined factor case. 
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3 

Y 

3.4 sample AH ALY SBS 

Reference (8) presents some results of the strength in- 
vestigation. Those results showed good agreement with pre- 
vious analyses in the literature. The present version of 
ST3EN was used to generate a failure locus of a 0/90 gra- 
phite-epoxy laminate at a temperature of 200 F below the 
curing temperature; for all stress states of the fora 
(Nx,Ny,Sxy,Hx,My,i1xy) = (Nx, Ny,0,0 /0,0,0,) . Since there are 
only two non zero stress components (Nx and Ny ) , the locus 
of failure stress states can be shown as a curve in Sx-Ny 
space. This curve is shown in Figure 15. This failure lo- 
cus was generated by running STREN for 37 combinations of Mx 
and Ny ( i.e., (Nx,>Iy) = (1.,0.) , (1 . , . 5) , . . . , etc. ) . 



= a + b - 



= - 1 . 
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Figure 27. Failure Locus for 0/90 Graphite-Epoxy Laminate 
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